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Abstract. In a recent publication [Y.D. Wang and A. A. Clerk, Phys. Rev. Lett. 108, 
153603 (2012)], we demonstrated that one can use interference to significantly increase 
the fidelity of state transfer between two electromagnetic cavities coupled to a common 
mechanical resonator over a naive sequential-transfer scheme based on two swap 
operations. This involved making use of a delocalized electromagnetic mode which 
is decoupled from the mechanical resonator, a so-called "mechanically- dark" mode. 
Here, we demonstrate the existence of a new "hybrid" state transfer scheme which 
incorporates the best elements of the dark-mode scheme (protection against mechanical 
dissipation) and the double-swap scheme (fast operation time). Importantly, this 
new scheme also does not require the mechanical resonator to be prepared initially 
in its ground state. We also provide additional details on the previously- described 
interference-enhanced transfer schemes, and provide an enhanced discussion of how 
the interference physics here is intimately related to the optomechanical analogue of 
electromagnetically-induced transparency (EIT). We also compare the various transfer 
schemes over a wide range of relevant experimental parameters, producing a "phase 
diagram" showing the the optimal transfer scheme for different points in parameter 
space. 
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1. Introduction 

The fields of optomechanics and electromechanics are interested in tlie pliysics of 
systems where a mechanical resonator is coupled to either a driven optical cavity, 
or a driven electronic circuit. The interaction of an optical cavity and mechanical 
oscillation due to radiation pressure or photon thermal effect has been studied in early 
experiments (see. e.g. [1,2]). In the past few years, these fields have achieved significant 
progress- ground state cooling and many-photon strong coupling have been realized in 
experiments on both optomechanical and electromechanical systems [3, 4, 5, 6, 7]. As 
mechanical resonators can achieve strong coupling with both microwave and optical 
cavities, this recent experimental progress may soon enable a powerful new application 
of these systems: the ability to transfer quantum states or traveling photons between 
optical cavities and microwave-frequency electrical systems. Such optomechanically- 
enabled state transfer could allow one to combine the complimentary advantages of both 
microwaves and optics. For example, one could imagine using a superconducting qubit 
to prepare a highly non-classical state of a microwave cavity (as recently demonstrated 
in Ref. [8]), and then transfer this state to an optical cavity, which provides an interface 
for long-distance information propagation via photons. Optomechanical state transfer 
could also serve as a key component in quantum information processing networks which 
combine a variety of different elements over long distances (see, e.g., [9]). 

In such schemes, one would be interested in the ability to transfer intra-cavity 
states, as well as the states of itinerant photons incident on one cavity. To be explicit, 
by intra-cavity state transfer we mean that the state to be transferred is the initial 
(t — 0) state of the first cavity; this could even be a non-classical state prepared by, e.g., 
having cavity 1 interacting with a qubit at t < 0. The goal is now to design a protocol 
so that cavity 2 ends by being in this same state. In contrast, for itinerant photon 
transfer schemes, the goal is to have a wave packet incident on cavity 1 be faithfully 
reproduced at the output of cavity 2. The bandwidth of this wave packet could be 
much much smaller than the damping rates of the cavities. Such conversion of itinerant 
photons of completely different wave length has various intriguing applications and has 
been discussed in several recent works [15, 17, 18]. In this paper, we will analyze both 
of the two transfer tasks. 

The problem of transferring quantum states between a mechanical resonator and 
cavity photons was the subject of several early studies [10, 11], utilizing strategies well- 
known in the field of cavity QED [12]. While the systems considered there involve 
only two resonators (an optical cavity and a mechanical resonator), similar approaches 
could be used in a three-resonator system where, e.g., both a microwave and an optical 
cavity are jointly coupled to a mechanical resonator; this is the subject of more recent 
studies [13, 14, 15]. For the problem of transferring intra-cavity states, these works 
consider the use of two successive "swap" operations in a two-cavity optomechanical 
system (Fig. 1). One pulses the optomechanical interactions such that initially, only 
the first cavity and the mechanical resonator are coupled. The coupling is then left 
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on for just long enough to exchange the states of the first cavity and the mechanical 
resonator. This is then repeated to exchange the mechanical and the second cavity 
states (i.e. the cavity 1 - mechanical resonator coupling is turned off, the cavity 2 - 
mechanical resonator coupling is turned on) [13, 14]. Note that single-swap operations 
between an optical cavity and mechanical resonator have recently been demonstrated 
experimentally [5, 16]. While straightforward, the "double-swap" scheme for intra- 
cavity state transfer has the drawback that as the state to be transferred occupies the 
mechanical resonator during the transfer protocol, it will be strongly degraded by any 
thermal noise driving the mechanical resonator. This remains true even if one first 
prepares the mechanical resonator in the ground state. 

In a recent work [17], we addressed the above problem by now looking at the 
two-cavity optomechanical system as a true three-mode system, as opposed to only 
considering two modes (i.e. one cavity mode and the mechanical resonance) at a time. 
One finds that there exists in general a delocalized cavity eigenmode of the coherent 
optomechanical Hamiltonian which is decoupled from the mechanical resonator, the so- 
called "mechanically-dark" mode. Using this mode enables high fidelity state transfer 
even in the presence of strong mechanical dissipation; similar ideas were also studied 
in [18]. In this paper, we provide additional technical details of the transfer protocols 
discussed in Ref. [17], in particular details on the analytic calculation of transfer 
fidelities, as well as a thorough discussion of the dependence of the transfer fidelity 
on the initial state of the mechanical resonator. Further, we introduce a new scheme to 
deal with one of the key difficulties of the adiabatic transfer scheme discussed in Ref. [17] , 
namely that to adiabatically use the dark mode for state transfer, the transfer protocol 
must be very slow, and hence is extremely susceptible to degradation from cavity losses. 
The new scheme we describe is a so-called "hybrid" scheme, as it only partially involves 
the dark-mode described above. While this means that it is not completely immune 
to mechanical heating effects, we show that it allows for a much faster transfer time, 
and is hence less susceptible to cavity losses. The hybrid scheme also has the surprising 
advantage that it does not require the mechanical resonator to be initially prepared in 
its ground state; we provide a heuristic explanation for this insensitivity in the text. As 
we analyze in detail in the paper, the hybrid scheme works best for typical experimental 
parameters where both cavity dissipation and mechanical noise are equally relevant. We 
give a comparison to the different features of each intra cavity state transfer scheme. 
We also discuss the transfer of itinerant photons and the underlying optomechanical 
EIT mechanism. 

The remainder of this paper is organized as follows. The basic system Hamiltonian 
and the definition of the dark mode are discussed in Sec. 2. In Sec. 3, we introduce and 
describe three basic schemes to implement intra-cavity state transfer: the double-swap 
scheme of Ref. [13] (which does not use the dark mode), the adiabatic passage scheme 
of Refs. [17, 18] (which entirely uses the dark mode), and the new hybrid scheme (which 
only partially uses the dark mode). The fidelity and features of each scheme is analyzed. 
We also compare the three schemes and provide a "phase diagram" showing the optimal 



Using dark modes for high-fidelity optomechanical quantum state transfer 



4 



transfer scheme for various choices of experimental parameters (optomechanical coupling 
dissipation, temperature, etc.). In Sec. 4, we change focus to now consider the state 
transfer of itinerant photons, that is wave-packets incident on one of the two cavities; this 
expands upon the analysis given in Ref. [17]. Particular attention is given to clarifying 
the connection between the high fidelities possible here and optomechanical analogue of 
electromagnetically-induced transparency (EIT) [19, 20, 21]. Conclusions and remarks 
are provided in Sec. 5. 

2. Basics of the two-cavity optomechanical setup 

2.1. Setup and Hamiltonian 

For concreteness, we consider an optomechanical system where a single mechanical 
resonator is simultaneously coupled to both an optical cavity and a microwave cavity via 
dispersive couplings (see Fig. 1(a)); particular experimental realizations are discussed in 
Ref. [14, 15]. The cavities interact with the mechanical resonator through the standard 
radiation pressure coupling (/i = 1) 

^inM ^ gMk [a + a^) , (1) 

where bi is the annihilation operator of cavity i {i = 1,2), a is the annihilation 
operator of the mechanical resonator, gi is the single-photon optomechanical coupling 
strength to cavity i, and the mechanical frequency is ujm- Note that to date, almost 
all experimentally studied devices are firmly in the limit of a weak single-photon 
optomechanical coupling, gi <^ the exception being experiments where the collective 
motion of cold-atoms in a cavity acts as the effective mechanical degree of freedom [22]. 

We further specialize to the standard situation where each cavity is strongly driven, 
resulting in a large average photon number which serves to enhance the effective cavity- 
mechanical resonator coupling. It is useful to separate the cavity state into a classical 
amplitude bi plus a deviation di which accounts for the effects of noise and the interaction 
with the mechanics: 

bi = e-'^^'^' (hi + d^) , (2) 

where ujd,i is the drive frequency of cavity z. The classical amplitude bi is simply 
proportional to the amplitude of the laser drive applied to cavity i; without loss of 
generality, we take it to be real and positive, di can be interpreted as the cavity 
i annihilation operator in a frame displaced by the classical amplitude. For a large 
drive bi ^ 1 and small single-photon optomechanical couplings, we can safely neglect 
terms in the displaced frame which are quadratic in the displaced cavity operators, 
resulting in purely linear coupling between each cavity and the mechanical resonator 
(see, e.g., [23, 24]). Finally, we work in an interaction picture with respect to the two 
cavity drives. The resulting Hamiltonian in this displaced interaction picture (in the 
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absence of any dissipation) reads: 

Hq = umcl^cl — I^A^rfJrf^ — Gi i^^^di + d\Gi^ (3) 



=1,2 



Here, = uja^i — Vti is the detuning of the drive appHed to cavity i and Vti is the 
cavity resonance frequency; we have assumed detunings near the red-detuned mechanical 
sideband of each cavity (i.e. A^ ^ —^m) and anticipated taking the good cavity hmit 
f^i ^ {i^i is the energy damping rate of cavity z), which allows us to make a 
rotating-wave approximation in writing the optomechanical interaction. The driven (or 
many-photon) optomechanical coupling between the mechanical resonator and cavity i 
is denoted as Gi = hiQi (we let > throughout the paper) ; note that hi is proportional 
to the drive amplitude applied to cavity and thus can be controlled in time. 

We see quite clearly from Eq. (3) that in our displaced frame, the optomechanical 
Hamiltonian has exactly the form we need for state transfer: the optomechanical 
interaction can move quanta from each cavity to the mechanical resonator, and vice- 
versa. We will thus be interested in quantum state transfer in the displaced frame of 
each cavity. Throughout this paper, the quantum state to be transferred is represented 
by di. This state of interest thus sits atop a large classical coherent state \hi) which 
is used to generate a large effective optomechanical coupling. Equivalently put, if we 
want to transfer a given state from cavity 1 to cavity 2, cavity 1 should be initially 
prepared in a state which is just the state displaced an amount hi in phase 
space. The role of this phase-space displacement is solely to achieve a large effective 
optomechanical coupling. Note that these two parts of the cavity state will be spectrally 
separated: the background coherent-state part of the cavity field used to generate a 
coupling is concentrated at the red-detuned sideband (i.e. Vti — oum)^ whereas the state 
to be transferred is located roughly within a cavity bandwidth k of the cavity resonance 
frequency Qi. 

We now turn to including dissipation in our system. The two cavities and 
mechanical resonator are assumed to be coupled to independent Ohmic baths that 
describe both internal loss and the couplings to the extra-cavity modes used to drive 
each cavity. For simplicity, we take each cavity to be single-sided (i.e. each cavity has 
a single input-output port). The master equation for the reduced density matrix p{t) 
describing both cavities and the mechanical resonator takes the form: 

y5 (t) = -I p (t)] + CcP {t) + Cup {t) (4) 

where the dissipation can be described by Lindblad type dissipator 

1=1,2 1=1,2 

CmP {t) = ^{Nm + 1)V (a) p{t)+l {Nm) V {a)) p (t) . (6) 

Here 7 (/^^) denote the energy decay rate of the mechanical resonator (cavity i), and 
Nm {Ni) denote the bath temperature, expressed as the number of thermal quanta, of 
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Figure 1. (Color online) (a) The two-cavity optomechanical system where a 
mechanical resonator (red bar in the middle) is coupled simultaneously to an optical 
cavity (the Fabry-Perot cavity in blue) and a microwave cavity (the LC circuit in 
green). Possible experimental realizations of this setup has been described in detail 
in Ref. [14, 15]. (b) The energy diagram of the three eigenmodes of the system for 
Gi = (left) and Gi ^ (right). Gi is the driven coupling strength of the mechanical 
resonator to the two cavities which is tunable with the drives (see Eq. (3) and the 
following text for more details). The middle red line corresponds to the mechanically- 
dark mode, whose frequency is independent of the coupling. ^/G1 -\- G2 is the energy 
difference between the dark mode and the two mixed modes. 



the bath coupled to the mechanical resonator (cavity i). The super-operator V{A) is 
defined as: 

V (i) p (t) = Ap (t) it - - jiti^ ^ (^) 1^ . (7) 

Finally, note that we are assuming throughout the paper that there is negligible phase 
noise in the large driving amplitudes on each cavity used to generate the optomechanical 
couplings Gi. The effect of such noise on coherent transfer was studied in Ref. [25, 26]; 
the resulting requirements are similar to those needed for cavity cooling [26], and are 
thus within reach of experiment. 



2.2. Mechanically- dark mode 

To understand how interference could be used to enhance mechanically-mediated state 
transfer in our system, we first note that the mechanical resonator only couples to 
a particular linear combination of the cavity modes. Defining new canonical mode 
operators 

and 

the Hamiltonian Eq. (3) can be rewritten as 



Using dark modes for high-fidelity optomechanical quantum state transfer 



7 



where we have assumed the two cavities are each driven at the red-detuned sideband 
= A = —cjm' We see that the "mechanicaUy-bright" mode described by Cbr is 
coupled with the mechanics, whereas the "mechanicaUy-dark" mode described by Cdk is 
decoupled from the mechanics. 

The above Hamiltonian is trivially diagonalized as 

i=+ ,-,dk 

where ^ 

2±W = ^(4rW±a), (12) 

describes hybridized modes formed from the mechanical resonator a, and the 
mechanically-bright cavity mode Cbr (t), and the mode eigenfrequencies are given by: 

uj^^ujM± ^JGl + Gl (13) 
^dk = (14) 

The energy structure of the three modes is shown in Fig. 1(b). 

As the dark mode Cdk(^) of the two cavities does not couple to the mechanical 
resonator (c.f. Eq. (10)), one would expect it to have immunity against the effects of 
mechanical dissipation (c.f. Eq. (6)). This simple fact will provide the foundation for 
several state transfer schemes that are robust against mechanical dissipation. However, 
it is worth noting that this protection is not exact when one considers the case of 
asymmetric cavity dissipation, that is K\ ^ One can easily understand this by 
just looking at the damping effects of the cavity dissipation. The Heisenberg-Langevin 
equation of motion for each cavity mode operator will have a damping term of the form: 

^rf. = (-zcjm - ^^2/2) di + ... (15) 

Some simple algebra now shows that if we re-write these expressions in the eigenmode 
basis, then for the damping terms in the equation of motion will mix the dark 

cavity mode and the bright cavity mode. The result is that mechanical noise can now 
corrupt the dark mode: mechanical noise will heat the coupled bright mode, and this 
noise will then drive the cavity dark mode via the cavity damping terms. We can see 
this more explicitly by writing the super-operator of cavity dissipation in the eigenmode 
basis: 

+ ^ ^ (^Cbr, C^k j P it) 

+ 2 r^k> Cbr j P [t) (16) 
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with the super-operator is defined as 



^ (i, b) p (t) ^Ap{t)B-- [bA, p it)] 



(17) 



The last two fines of Eq. (16) represent the mixing of the bright and dark mode; these 
terms vanish when both ni = and A^i = A^2- This mixing effect leads to an effective 
heating of the dark mode which is a factor ^ {tii — f^2)/G smaller than the heating of the 
mechanically-bright cavity mode. While the mechanically-dark mode is only completely 
immune from mechanical noise for completely symmetric dissipation, we will show that 
schemes utilizing it can still be highly effective even when there are deviations from this 
perfect symmetry condition. Note that a definition of dark mode for Ki ^ K2 based on 
a perturbative treatment of the cavity damping asymmetry was discussed in [18]. 

3. Intra- cavity state transfer 

We first consider protocols which allow the two-cavity optomechanical system in Fig. 1 
to be used to transfer intra-cavity states. The initial state of interest is prepared in 
cavity 1, and the goal is to have this state end up in cavity 2 at the end of the transfer 
protocol with a maximal fidelity. We will discuss 3 different schemes to implement 
such a transfer process: the double-swap protocol (which does not use the dark mode), 
the adiabatic passage protocol (entirely based on the dark mode) and the hybrid swap 
protocol (which only partially uses the dark mode). In the absence of any cavity or 
mechanical dissipation, all these schemes allow transfer with perfect fidelity. However, 
in the presence of dissipation, this is not so; moreover, each scheme has its own relative 
merits and disadvantages versus the different kinds of system dissipation. We will discuss 
this in what follows, with the goal of providing insight into what the optimal strategy 
is given a certain set of system parameters. 



In order to quantify the performance of a transfer protocol, we consider the 
Ulhmann fidelity F [27]. 



where pi (pf) denote the density matrix of initial state to be transferred (the final state 
in cavity 2). Note that we will optimize the fidelity over simple rotations in phase space 
(so that if pf is a rotated version of pi, F = 1). 

We will begin by considering the simple case of transferring a pure Gaussian state 
between the two cavities (i.e. a pure state whose Wigner function is Gaussian). For 
such states, simple analytic expressions can be derived for the transfer fidelity [17, 28], 
allowing quantitative insight into the various physical processes which degrade transfer. 
Starting with an initial pure Gaussian state in cavity 1, the transfer fidelity F takes the 
general form (for details, see Appendix A): 




(18) 




(19) 
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Figure 2. (Color online) Schematics for the three intra-cavity state transfer schemes 
considered in this paper. Top panel: the curves in the boxes show the required time- 
dependences of the two optomechanical couplings in each scheme: (a) double-swap 
(b) adiabatic passage, and (c) hybrid scheme. Gi {i = 1,2) is the optomechanical 
coupling to cavity i. The brown dashed line denotes the maximum value of the 
coupling strength Gi^max- Lower panel: schematics illustrating how the transfer works 
in each scheme. Here the blue, red and green dots represent (respectively) cavity 1, the 
mechanical resonator, and cavity 2; the wave packet in cavity 1 represents the state 
to be transferred. Arrows indicate the optomechanical couplings, with the magnitude 
indicating their magnitude. 

F depends on just two positive-definite parameters nh and A. nh can be regarded as 
an effective number of thermal quanta which quantifies the heating of the state during 
the transfer protocol. This heating can arise both from noise emanating from the cavity 
and mechanical dissipative baths, as well as from any initial thermal population in the 
mechanical resonator. In contrast A characterizes the deleterious effects of amplitude 
decay during the state transfer protocol, i.e. the decay of the average of {d). Efficient 
transfer requires minimizing both effects. 

The values of nh and A depend both on the nature of the initial Gaussian state, and 
on the precise nature of the time evolution associated with the transfer protocol. The 
3 possible protocols correspond to different time-dependences of the optomechanical 
couplings Gi{t) and ^2(^)5 as shown in Fig. 2. In the following subsections, we will 
discuss the behavior of fih, A and the fidelity F for each of the three intra-cavity transfer 
schemes. 

3.1. Double-swap protocol 

As already discussed in the introduction, the most straightforward manner to exploit 
Eq. (3) for state transfer is via two sequential swap operations (see Fig. 2(a)). As this 
scheme only involves at most two modes interacting simultaneously, the mechanically- 
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dark mode introduce earlier plays no role in its physics. In order to provide a basis 
for comparison for the dark-mode schemes that follow, we derive both the fidelity for 
this scheme, as well as discuss why this scheme requires the mechanical resonator to be 
initially prepared in its ground state. 

Using the Hamiltonian in Eq. (3), the double-swap protocol requires that we 
first only turn on the interaction between cavity 1 and the mechanical resonator (i.e. 
Gi (t) = Gi,max and G2 {t) = as shown in Fig. 2(a)). In the absence of dissipation, the 
time evolution of cavity 1 and mechanical mode reads 

a (t) = e-'^^' [a (0) cos(Gi,^axi) - idi (0) sin(Gi,^axt)) 

d, (t) = e-^-^^ (di (0) cos(Gi,„,axt) - id (0) sin(Gi,^axt)) (20) 

It thus follows that at a time tis — 7r/(2Gi^niax), cavity 1 and the mechanical resonator 
will have perfectly exchanged their states (up to a trivial phase factors) [10]: 

a ihs) = e-^^(*-)cli (0) , d, ihs) = e-^^(*-)a (0) (21) 

with 9{t) = ujMt + 7r/2. As the initial cavity-1 state is now in the mechanical resonator, 
one next turns off Gi and turns on G2 = G2,max for a time = 7r/2G2,max to similarly 
exchange the mechanical and cavity-2 states (see Fig. 2(a)). At the end of the operation, 

di{h, + hs) = e-'^(*^=+*^=)J2(0), 4(tis + h,) = e-^^(*i=+*^=) Ji(0) (22) 

Thus, without dissipation, one can perfectly transfer an arbitrary cavity-1 state to 
cavity-2. 

In the presence of dissipation, one might expect that high fidelity transfer is still 
possible if the total transfer time ti^ + is much shorter than the relevant dissipative 
time scales. This expectation can be easily quantified in the case of Gaussian state 
transfer, using the general expression for the transfer fidelity F given in Eq. (19). 
Not surprisingly, one finds that each cavity must achieve the many-photon strong 
coupling condition i^i <C Gi,max; otherwise, the decay of the state in the two cavities 
during the transfer makes high fidelity impossible. In this strong coupling regime, the 
heating parameter fih and amplitude-decay parameter A that determine the fidelity take 
particularly simple forms. The full expression for an arbitrary squeezed state (along 
with calculation details) is presented in Appendix B.l. Here, we present the results for 
the strong coupling case with a coherent state \a) as the initial state. We also assume 
the mechanical resonator has initially been prepared in the ground state. One finds that 
to leading order in (k + 7)/G: 

^ Us ^ 7^M^, 

i 

^ = HE^T^^-^l"l2 2G- ^^^^ 

i 

The effective number of thermal quanta nh generated during each time interval 
tis is simply given by the duration of the interval times the average heating rate 




Figure 3. (Color online) The performance of the double-swap scheme versus coupling 
strength G = Gi,max = G^2,max, for the transfer of a coherent state having amplitude 
a = 1. Shown are (a) transfer fidelity, (b) heating parameter nh (note that a 
logarithmic scale is used), (c) damping parameter A. nh and A are defined after 
Eq. (19). We have fixed the bath temperature to be T = 1.5 K, (which would 
be compatible with the use of a superconducting microwave cavity fabricated from 
niobium), and taken = {hzi -\- hZ2) /"^ = 27r x 50 KHz (which is appropriate to the 
system envisaged in Ref. [14]). The blue solid lines correspond to hzi = and a 
mechanical resonator which is initially prepared in its ground state; the red dashed 
lines correspond to hzi = K2 but without any mechanical pre-cooling (i.e. the mechanical 
resonator is initially in thermal equilibrium at the bath temperature). The green 
dashed-dotted line corresponds to the case where the mechanical resonator is precooled, 
but i<i2 = 4/^1. Note that all lines coincide in plot (c), as the amplitude decay parameter 
A is insensitive to asymmetry in the cavity damping, and to mechanical pre-cooling. 
The other parameters are 7 = 27r x 1 KHz, ujm = — = 27r x 10 MHz. Cavity 1 (2) 
is a microwave (optical) cavity: Cli/27T = 10 GHz (1^2 /27r = 100 THz). The transfer 
time is taken to be the optimal time 4 as defined in the text. 



{jNm + i^iNi) /2. Similarly, the amplitude decay A in each interval is just the duration 
of the interval, times the amplitude of the initial state, times the average amplitude 
decay rate (/^^ + 7)/4. The last approximation in each equation above corresponds to 
the limit of symmetric cavity parameters /^^ = 7, = G and A^^ 0. 

The above result implies that for mechanical heating to not be a problem, one 
needs that the optomechanical couplings G not only be larger than but also larger 
than the mechanical thermal decoherence rate ^Nm- This regime was recently realized 
in experiment [5]. Plotted in Fig. 3 are transfer fidelities corresponding to the transfer 
of a coherent state between a microwave cavity and an optical cavity, using numbers 
appropriate to kind of setup envisaged in Ref. [14]; the behavior of the parameters nh 
and A are also shown over a range of parameters. We have taken a relatively high bath 
temperature of 1.5 Kelvin, as from an experimental point of view, it would be highly 
advantageous to not have to cool the optical cavity to tens of mK. We note that such a 
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high operating temperature is compatible with the use of a superconducting microwave 
resonator fabricated from niobium. For such a high temperature, one clearly sees that 
the double-swap scheme is only able to yield good fidelities at extremely strong values 
of the optomechanical coupling G. 

Degradation due to mechanical heating is even more severe when considering the 
transfer of squeezed states (c.f. Eq. (B.16)): this enhanced sensitivity is due to the state 
having narrower features in phase space. As shown in Eq. (B.18), when transferring a 
squeezed vacuum state, the general expression of Eq. (19) for the fidelity still holds, but 
the heating factor fih is enhanced exponentially compared to the coherent-state case. 
This is a direct consequence of the initial state having extremely narrow features in phase 
space; it suggests that the double-swap scheme will be especially poor in transferring 
such states (in comparison to schemes that attempt to suppress mechanical heating). A 
similar conclusion can be expected when considering the transfer of non-classical states 
(e.g. a Fock state), which also have fine structures in phase space. 

We thus see that while the double-swap protocol works well at extremely low 
temperatures, it is already problematic at moderate values of mechanical bath 
temperature. Given the experimental advantages of being able to use such temperatures, 
there is ample motivation to consider alternate transfer protocols which are more robust 
against the effects of mechanical thermal noise; we do this in the following subsections. 

Before proceeding, we consider the importance of precooling to the double-swap 
scheme. In Eq. (23), we have assumed the optimal situation where the mechanical 
resonator is initially in its ground state [13, 14, 17], which could be achieved by initially 
swapping a cavity ground state into the mechanics, or by using conventional cavity 
cooling [23] . If in contrast one does not initially cool the mechanics to the ground state, 
its initial thermal population Nm,o will make a significant additional contribution to the 
heating parameter nh- One finds: 

= 2^ ^ Us + ( — 1 Nm,o- (24) 

One might naively think that the initial thermal population of the mechanical resonator 
should be irrelevant, as it will just be swapped to cavity- 1 in the initial swap operation 
of the protocol. However, this is not the case due to cavity dissipation. We can obtain a 
simple understanding of this by making an analogy between our system and the evolution 
of a quantum two-level system (TLS) in a magnetic field. Including only the damping 
effects of dissipation, and letting v — ^a(t), Ji(t)^ , the Heisenberg equation of motion 
for the cavity- 1 and mechanical mode operators during the first swap pulse takes the 
form 

J^v^-^UMV-^v-^G,,^..^^ j " " ^ -1 ^''^ 

If we now interpret i; as a spinor describing a quantum two-level system written 
in a basis of eigenstates, then for Ki =7 = 0, the above equation corresponds to 
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Figure 4. (Color online) (a) Bloch sphere representation of a two-cavity swap 
operation. Each swap operation in the double-swap transfer scheme (c.f. Sec. 3.1) can 
be understood in analogy to the evolution of a two- level system in a magnetic field; 
one takes the mechanical mode operator d and cavity mode operator di to correspond 
to the az eigenstates of the TLS. Without dissipation, each swap operation then 
corresponds to a tt rotation (as indicated by the blue circle on the sphere) around the 
effective magnetic field Bq^ which is oriented along the x axis and is proportional to the 
optomechanical coupling Gi^max- As discussed in the text, including dissipation changes 
this picture. The large asymmetry between the cavity damping Hi and mechanical 
damping 7 tilts the effective magnetic field towards z-direction to B'^. A rotation 
about this tilted field now cannot perfectly rotate d into di (as shown by the brown 
circle on the sphere): a perfect swap is thus not possible. This explains why the the 
double-swap scheme is so susceptible to an initial mechanical thermal population, (b) 
A similar Bloch sphere representation can be used to understand the hybrid transfer 
scheme of Sec. 3.3; now however the eigenstates of our effective TLS correspond to 
the mechanical mode and the "bright" cavity mode Cbr (c.f. Eq. (8)), and the effective 
rotation that is needed is a 27r rotation (see discussion after Eq. (29)). As such a 
rotation can be performed even with a magnetic field that is not in the equatorial 
plane, the hybrid scheme is insensitive to the initial mechanical population. 



the precession of the TLS in a magnetic field of strength oc Gi^max oriented in the x 
direction. The perfect swap of cavity and mechanical states occurring at a time t = tis 
(c.f. Eq. (21)) thus corresponds directly to a tt rotation of the effective TLS in the Bloch 
sphere, as sketched in Fig. 4(a) with the blue circle. Turning now to the dissipative terms 
in Eq. (25), we see that there are two effects. The first is an overall exponential decay 
induced by the average damping rate (/^i + 7)/2; this corresponds to the second term 
on the RHS of Eq. (25). This decay on its own would still allow for a perfect swap, 
in that a (tig) would only depend on rfi(O) and not a(0). In contrast, the last term on 
the RHS of Eq. (25) (proportional to (/^i — 7)/2) is like an effective z-magnetic field 
acting on the TLS (albeit with an imaginary magnitude). This terms makes a perfect 
swap impossible, as the effective magnetic field acting on our TLS is no lo nger purely 
transverse, i.e., in Fig. 4(a), the rotation axis is changed from to S^. It is this 
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Figure 5. (Color online) The performance of the adiabatic passage scheme versus the 
maximum coupling strength G Gi max 

Go 

max foi" the transfer of a coherent state 
|a = 1), using parameters is identical to Fig. 3. Shown are: (a) Fidelity, (b) the heating 
parameter nh (notice that a logarithmic scale is used in this plot), (c) the damping 
parameter A. The blue solid lines correspond to tzi = K2 and a mechanical resonator 
that is initially in the ground state; the red dashed line corresponds to the same without 
precooling. The greed dashed-dotted lines correspond to a pre-cooled mechanical 
resonator but unequal cavity damping K2 = 4/^i and hzi -\- tZ2 = 2/^. Even though 
asymmetrical cavity damping destroys the perfect protection of the mechanically-dark 
mode from mechanical dissipation, it only results in a small decrease of the fidelity of 
the adiabatic transfer protocol; see text for details. For each value of the coupling, we 
have optimized the speed /3 and final time tf of the protocol to obtain a maximum 
fidelity (see main text). 

term that makes the double-swap protocol extremely susceptible to any initial thermal 
population in the mechanics: this population is not fully swapped away in the first 
swap operation, and thus ends up being transferred to the final state of cavity 2. As 
shown in Fig. 3 by the dashed line, if one does not precool, the transfer fidelity of the 
double-swap scheme is greatly suppressed. As we will see in Sec. 3.3, there does exist 
a transfer scheme (the hybrid scheme) which does not require any initial precooling of 
the mechanical resonator. 

3.2. Adiabatic passage 

The double-swap scheme of the previous section only involves two modes interacting 
with one another at any given time, and hence could not take advantage of the full 
structure of our three-mode (two cavities, one mechanical resonator) system. Here, we 
consider an alternative scheme (first discussed in Refs. [17, 18]) which explicitly makes 
use of the "mechanically-dark" mode Cdk introduced in Sec. 2.2 (c.f. Eq. (9)), a mode 
which is not coupled directly to the mechanical resonator and hence protected against 
mechanical dissipation. The basic idea is to use this dark mode to perform the state 
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transfer from cavity 1 to 2 entirely via the dark mode. This is done by adiabaticaUy 
modulating the optomechanical couplings Gi, G2 in time, so that at t = the dark mode 
is equal to the cavity- 1 mode, and at the end of the protocol t = it is equal to the 
cavity-2 mode. This involves using a pulse sequence where Gi and G2 are modulated 
in time as shown in Fig. 2(b). Gi{t) is increased from to Gi^max, while G2 (^) is 
simultaneously decreased from G2,max to 0. For the transfer to be perfectly adiabatic, 
one requires that the transfer speed be slow compared to the energetic gap separating 
the dark mode and the two coupled modes c± defined in Eq. (12). This energetic gap is 
given by ^yGl + G^ (c.f. Eq. (13)). If one is indeed adiabatic, the dark mode will evolve 
during this protocol from being —di at t = to ^2 at the end of the protocol at a time 



The adiabatic passage scheme described here is analogous to the well known 
Stimulated Raman Adiabatic Passage (STIRAP) scheme employed in atomic physics 
(see Ref. [29] and reference therein). In STIRAP schemes, one employs a dark-state of 
a 3-level atom, and also uses a similar "counter-intuitive" pulse sequence to transfer an 
atomic population between two states. In contrast, here we have an entire manifold of 
dark states (i.e. any states in the subspace generated by c^j^), and thus the potential to 
transfer an arbitrary cavity state using the adiabatic protocol. 

To quantify the effectiveness and fidelity of this scheme, we again consider the 
transfer for a Gaussian state; the fidelity F thus takes the general form described by 
Eq. (19). While there are many possible ways to implement the adiabatic passage 
scheme, for simplicity, we use throughout the simple pulse shapes 



This form keeps the energy splitting between the dark and coupled modes constant (i.e. 
Gl (t) + Gl (t) = G^), a feature that has been argued to be optimal [30]. The parameter 
/? represents the overall speed of the modulation; one would require /3 <C G to be in the 
purely adiabatic limit. 

Shown in Fig. 5 are results for fidelity, heating parameter fih and amplitude decay 
parameter A for the adiabatic passage scheme, using parameters identical to Fig. 3 (in 
particular, a bath temperature of 1.5 K). For each data point, we numerically found the 
optimal values of the transfer speed /3 and transfer time t/, and use these to obtain 
the fidelity. As expected, using the dark-mode leads to a dramatically suppressed 
nh compared to the double swap scheme (i.e. compare solid curves in Fig. 3(b) and 
Fig. 5(b)). If the mechanical bath temperature is sufficiently high, the adiabatic transfer 
scheme yields far higher fidelities than the double-swap scheme. This behavior is shown 
explicitly in Fig. 6, where the transfer fidelity versus temperature for both these schemes 
are plotted. We will quantify this advantage in more detail in Sec. 3.4, where we compare 
and contrast the performance of all transfer schemes. 

One might be puzzled by the results of Fig. 5(b), which shows that the effective 
heating parameter nh is never exactly zero, even in the case /^i = /^2, where the dark 
mode is completely isolated from the mechanical resonator. The reason is simple to 



t^tf. 




(26) 
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understand: because of the effects of cavity damping, the optimal pulse speed /3 in the 
adiabatic transfer protocol can never be zero, and hence for realistic cavity parameters, 
one will never be in the perfectly adiabatic limit. The result is some mixing into the non- 
adiabatic subspace, and hence some mechanically-induced heating. If one were in the 
adiabatic limit /3 <C G, mechanical noise would indeed not be a problem, but the slow 
transfer time would mean that the decay of the state in the two cavities would greatly 
degrade the fidelity. Thus, the two confiicting requirements of better adiabaticity (to 
suppress mechanical heating) and faster operation time (to suppress cavity damping 
effects) result in an optimal, non-zero value of /?opt and an optimal transfer time topt- 
High fidelity transfer is only possible when G ^ /3 ^ a^, as it is only in this regime where 
one can both be reasonably adiabatic as well as reasonably fast compared to the cavity 
damping rate One also finds that the optimal value /3opt has a marked dependence on 
temperature, as the relative importance of being adiabatic increases as the mechanical 
thermal noise increases. The fact that one is not in the perfectly adiabatic limit also 
leads to the requirement that the mechanical resonator to initially start in the ground 
state, as shown by the marked difference between the red-dashed and blue-solid lines in 
Fig. 5 (a) and (b). 

The above discussion also explains why for realistic parameters, there is no stringent 
requirement that the cavity damping be symmetric, i.e. — 1^2] even highly asymmetric 
damping only slightly reduces the transfer fidelity (see the green dash-dot line in Fig. 5). 
As discussed in Sec. 2.2, asymmetric cavity damping mixes the cavity dark and bright 
modes, thus opening up the dark mode to mechanical noise. However, as the optimal 
transfer speed never corresponds to being in the perfectly adiabatic limit, the additional 
exposure to the mechanical noise caused by unequal has rather small contribution to 
the heating parameter fih above that stemming from non- adiabatic transitions. The 
upshot is that the adiabatic transfer scheme can be useful even if one has rather 
dissimilar cavity dissipation rates. 

Finally, we note that in the perfectly adiabatic limit, the parameters nh and A can 
be calculated analytically [18]. However, for moderate G/k ratios, the optimal pulse 
shapes are not perfectly adiabatic, and hence these expressions are not particularly 
relevant. 

3.3. Hybrid transfer scheme 

So far in this section, we have considered two very different approaches to state transfer, 
each effective in a very different limit. The double-swap scheme of Sec. 3.1 does 
not involve the mechanically-dark mode in any way, and is hence more susceptible to 
mechanical thermal noise; however, it is relatively fast (the transfer time t oc 1/G), and 
hence less susceptible to degradation due to cavity damping. In contrast, the adiabatic 
passage scheme of Sec. 3.2 attempts to only use the dark mode for state transfer. The 
result is a strong resilience against mechanical noise. However, the constraint of being 
near-adiabatic means that this scheme is slow (i.e. the transfer time t ^ 1/G), and 
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Figure 6. (Color online) Comparison of the performance of the three intra-cavity 
transfer schemes for the transfer of a coherent state of amplitude a = 1, as a function 
of bath temperature. We have taken Gi^max = G^2,max = 27r x 2 MHz, the other 
parameters are identical to Figs. 3 (a). Shown are (a) fidelity, and (b) heating 
parameter nh- The blue dashed curve is the adiabatic transfer protocol, the green solid 
curve the hybrid scheme, and the red dash-dotted curve the double-swap scheme. As 
expected, the use of the mechanically-dark mode in the hybrid and adiabatic transfer 
schemes afford a much stronger resilience against the effects of mechanical heating. 

hence more susceptible (relative to the double-swap scheme) to the effects of cavity 
damping. 

Given that both schemes have their relative merits, one might be tempted to develop 
a strategy that incorporates the best features of both. We present such a scheme 
here, the so-called "hybrid" scheme, which utilizes equally the mechanically-dark mode 
(c.f. Eq. (9)) and mechanically-coupled cavity modes (c.f. Eq. (12)) for a state transfer 
between the two cavities. As we will see, the fact that half the transfer is done using 
the dark mode reduces the effects of mechanical thermal noise by a factor of two over 
the double-swap scheme. However, unlike the adiabatic passage scheme, this protocol is 
relatively fast (i.e. the transfer time is ^ l/^)' hence the effects of cavity damping 
are similar to what one would expect in the double-swap scheme. As an added bonus, we 
also find that this hybrid scheme does not require any precooling of the mechanics: high 
transfer fidelities are possible without having to initially cool the mechanical resonator 
to its ground state. For typical experimental situations where both cavity damping and 
mechanical heating are of equal importance, this hybrid scheme will in general yield 
the highest transfer fidelities (something we will make more clear in Sec. 3.4, where we 
compare the various intra-cavity state transfer protocols). 

The hybrid scheme we have in mind is realized in an extremely simple fashion: 
one simply turns on simultaneously both optomechanical couplings Gi and G2 to equal 
magnitudes, as shown in Fig. 2(c). It is simplest to consider the physics using the 
basis of dark and bright modes to describe the two cavities (c.f. Eqs. (8), (9)). Taking 
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<^i,max = <^2,max = G, One finds: 

(^1 = (cbr + Cdk) (27) 

4 = (cbr - Cdk) (28) 

The dark mode is uncoupled to the mechanical resonator, whereas the bright cavity 
mode and mechanical resonator evolve according to the Heisenberg equation of motion: 

where w = (cbr,S), and we have taken K.i = k and neglected noise terms for simphcity. 
Note the similarity between this equation and Eq. (25) which describes the first step of 
the double-swap scheme. 

Without dissipation (i.e. /^^ = 7 = 0), the situation here is analogous to the double- 
swap scheme: the evolution of the two interacting modes is equivalent to that of a TLS in 
a transverse magnetic field. In this case however, the cr^-eigenstates of the eflFective TLS 
corresponds to the cavity bright mode and the mechanical mode. From Eqs. (27), (28), 
it follows that to swap the states of the two cavities, we simply need to change the sign 
of the bright mode relative to the dark mode. This is easily achieved by performing a 
27r rotation of our effective TLS. Thus, after a time t = ths = ^r/ (a/2G), we find: 

rfi(ths) = e-''d2 (0) , d2(ths) = e-^'d, (0) (30) 

with 9 = ujMths + TT. Thus, similar to the double-swap scheme, the state of the two 
cavities are swapped. 

Unlike the double-swap process, this hybrid scheme involves both cavities being 
simultaneously coupled; this will have implications in terms of how the transfer is now 
affected by dissipation. On a heuristic level, the initial state in cavity 1 corresponds 
to having an equal initial population of the cavity bright and dark modes. Thus, 
half of the transfer is performed via the dark state, and we expect a better resilience 
against mechanical thermal noise than the double swap scheme. To quantify this, we 
again consider using the hybrid scheme to transfer a Gaussian state, including the full 
effects of cavity and mechanical dissipation. The fidelity F thus takes the general form 
given in Eq. (19) (see Appendix B.2 for details), with the heating and amplitude decay 
parameters in the strong coupling regime: 

— ^^hs 



4 2x/2 2G' 

A= H— ths-H^^. (31) 

For simplicity, we have taken the two cavities to be identical: Ki = ^2 = k^ 
Ni = N2 = Nc. The last approximation in each equation corresponds to the usual case 
where the mechanics dominates the heating, while the cavity dominates the amplitude 




Figure 7. (Color online) The performance of the hybrid scheme versus coupling 
strength G = Gi,max = G2,max for the transfer of a coherent state of amplitude a = 1, 
at a bath temperature T = 1.5 K; parameters are the same as Figs. 3 and 5. Shown 
are: (a) fidelity, (b) the heating parameter nh, (c) the damping parameter A. The blue 
solid line corresponds to /^i = /t:2 = and a mechanical resonator initially prepared 
in the ground state; the red dashed line corresponds to the same, but without any 
mechanical precooling (the mechanics starts from a thermal equilibrium state). As 
discussed in the main text, the hybrid scheme is almost completely insensitive to any 
initial thermal population in the mechanical resonator. The green dashed-dotted lines 
correspond to the case with pre-cooling but unequal cavity damping = 4/^i. The 
transfer time is taken to be the optimal time t^s as defined in the text. 



decay. Comparing the expression for nh with that of the the double swap case in 
Eq. (23), we find that the effect of mechanical thermal decoherence is reduced here by 
a factor of 2^/2. This comes from the fact that half of the input state is transferred via 
the dark mode, and from the shorter transfer time in the hybrid scheme, i.e., ti/{^/2G) 
versus tt/G. The effect of cavity damping is only slightly stronger than in the double 
swap case (by a factor 3a/2/4 ?^ 1.06). Fig. 6 shows the transfer fidelity of a a = 1 
coherent state versus temperature for both the hybrid and double-swap schemes, while 
Fig. 7 shows how the hybrid-scheme fidelity varies with coupling; both demonstrate the 
greater resilience against thermal heating of the hybrid scheme over the double-swap 
scheme. 

As already mentioned, an additional (and surprising) advantage of the hybrid 
scheme over both the adiabatic and double-swap schemes is that it does not require the 
mechanical resonator to start in the ground state. We can understand this heuristically 
by again using our analogy to the precession of a TLS. As we have explained, both 
the double-swap and the hybrid schemes corresponds to effective rotations in the Bloch 
sphere around a magnetic field Bg (see Fig. 4). In the absence of dissipation, this 
magnetic field is aligned along the x axis and can be used to perform a perfect tt rotation 
(in the double-swap scheme) or perfect 27r rotation (in the hybrid scheme). Including 
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mechanical and cavity dissipation induces new terms in the equations of motions, which 
can be interpreted as an additional (unwanted) magnetic field along the z-direction 
(c.f. Eqs. (25) and (29)). As discussed, this misalignment of the magnetic field means 
that a perfect tt pulse is no longer possible, and hence in the double-swap scheme, one 
is sensitive to the initial mechanical state. In contrast, in the hybrid scheme we need 
to perform a 27r rotation. This is possible even if the dissipation causes the effective 
magnetic field to be misaligned and not purely transverse. The result is that even with 
dissipation, the final cavity-2 state will remain completely independent of the initial 
mechanical thermal population. This is the reason that the performance with pre- 
cooling (blue solid line) and without pre-cooling (red dashed line) overlap as shown in 
Fig. 7. 

Finally, from Fig. 7, we also see that the hybrid scheme is insensitive to the 
asymmetry of cavity damping (the green dashed line only deviates from the blue solid 
line slightly). This is because the thermal heating here is dominated by the the use of 
the coupled cavity modes; the small heating of the dark mode that occurs when /^i ^ 
is only a small correction to this, especially given the relatively short transfer time of 
the scheme (it is the fastest of the three schemes we consider) . 

3.4' Comparison of intra- cavity state transfer protocols 

Having described the three different schemes for intra-cavity state transfer sketched 
in Fig. 2, we can now turn to compare them and discuss their relative merits. We 
stress that in the absence of any cavity or mechanical dissipation, all three schemes are 
able to achieve perfect state transfer of an arbitrary state. However, in the presence of 
dissipation, this will not be the case; moreover, the fidelity of each scheme will depend on 
the particulars of the dominant dissipation as well as magnitude of coupling strengths. 

To allow for an easy visual assessment of this competition, we have plotted in 
Fig. 8(a) and (b) "phase-diagrams" which indicate the optimal scheme for transferring 
di \a — 1) coherent state, as a function of both environmental temperature (taken to be 
the same for cavity and mechanical baths) and cavity damping rate These two axes 
characterize the dominant dissipation mechanics: heating from the mechanical bath, and 
amplitude decay due to the cavity damping. A different color (as indicated in the figure) 
is used to represent each of the three schemes; for each point, the colour corresponding 
to the highest-fidelity scheme is plotted, with an intensity proportional to the difference 
of the highest fidelity between the average fidelities of all 3 schemes. For simplicity, we 
have use symmetric cavity parameters Gi^max G2,max = G and Ki = K2 = 

The bottom left corner of both Fig. 8(a) and (b) corresponds to the limit where 
dissipative effects are minimal (i.e. low bath temperature and low cavity damping); as 
expected, all three schemes have almost equal performance. As we move away from the 
origin, dissipation is increased, and the near-equivalence of all three schemes in terms of 
fidelity is lost. As dissipation is increased, the fidelity of all the schemes considered will 
decrease, as shown in Fig. 8 (c) and (d). Note that all the schemes considered here take 
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Figure 8. (Color online) (a) A "phase diagram" showing the transfer protocol which 
yields the highest fidelity for the transfer of a coherent state of amplitude a = 1, as a 
function of environmental temperature and cavity damping, for a fixed optomechanical 
coupling Gi^max = G^2,max = 27r X 2 MHz, and for 7 = 27r x 100 Hz . Red points 
denote the double-swap protocol, green points the hybrid scheme, and blue points the 
adiabatic passage scheme. The saturation of each color is proportional to the fidelity 
difference between the optimal scheme and the remaining two schemes. We only show 
the points where the transfer fidelity F > 0.60. We also take ujm/'^'t^ = 10 MHz, 
(^i/27r = 10 GHz, and = 100 THz. (c) Contour plot showing fidelity of the 

optimal schemes, same parameters as (a), (b) and (d): Same as (a) and (c), except 
now the mechanical damping is increased to 7 = 27r x 1 KHz. 



a time at least on the order of 1/G to complete, and hence require the many-photon 
strong coupling condition G > tv to achieve a good fidelity. However, within this strong 
coupling regime, the performance of each scheme varies due to their different sensitivity 
to dissipation mechanisms: 

• The double-swap scheme (red dots in Fig. 8) is fast but it swaps the entire state 
into the mechanical mode as an intermediate step. Thus it is the scheme most 
fragile to the effects of mechanical heating, but the most robust against the effects 
of cavity damping. For an experimental setup that is able to achieve extremely 
low temperatures but cannot go deep into the strong coupling regime G ^ the 
double swap protocol will be the most suitable protocol to follow, as shown in the 
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upper left region of Fig. 8. Note from Eq. (23) that the importance of amphtude 
decay increases as one increases the average amphtude of the state to be transferred; 
this suggests the double-swap scheme would also have an advantage for the transfer 
of such states. On the other hand, the importance of mechanical heating becomes 
even stronger for squeezed states; we expect this to be also true for non-classical 
states having fine features in phase space. Thus, the double-swap scheme will be 
less advantageous for transferring such states. 

• In contrast, the adiabatic scheme (blue dots in Fig. 8) is largely immune to 
mechanical thermal noise; this protection is due to its use of the mechanically- 
dark mode. It is however easily degraded by cavity damping due to its necessarily 
slow transfer time ^/G). It is the best strategy when the effects of mechanical 
heating completely dominate the effects of cavity decay, as shown in the lower 
right region of Fig. 8. As discussed above, the heating effect becomes even more 
deleterious if transferring states with fine structure in phase space (such as squeezed 
states); the adiabatic scheme can thus be expected to be more optimal than the 
double-swap scheme for transferring such states. 

• The hybrid scheme (green dots in Fig. 8) utilizes both the mechanically-dark mode 
and the coupled-cavity modes for state transfer. The resulting fast transfer time 
and partial immunity to mechanical thermal noise make it the optimal scheme 
when both temperature and cavity damping are equally problematic, as shown in 
the intermediate region of Fig. 8. In addition, the hybrid scheme does not require 
one to initially prepare the mechanical resonator in the ground state; this could 
also present a strong practical advantage. 

4. Itinerant state transfer 

The previous section was entirely devoted to three schemes for transferring an intra- 
cavity state between two cavities. The state to be transferred is first prepared in one 
cavity, and then (via the interaction with the mechanical resonator) is transferred to the 
state of the second cavity. While this could have many useful applications (as already 
discussed), we have shown that all the intra-cavity transfer schemes require one to at 
least approach the many-photon strong-coupling condition, G > k.. In this section, 
we now switch gears and consider a state-transfer task that is possible even if one is 
not in this strong-coupling regime. This alternative task is the transfer of itinerant 
photons, i.e. converting an incoming wave-packet on cavity 1 into an outgoing wave- 
packet leaving cavity 2. Such conversion of itinerant photons is also something that 
could be of extreme utility, and has been discussed in several recent works [15, 17, 18]. 
Experimental demonstration based on mechanical resonator and a single optical cavity 
has been carried out recently [31]. We expand here on the discussion in [17], providing 
additional details and heuristic explanations for the phenomena which make itinerant 
state transfer possible in the two-cavity optomechanical system under consideration. 
We start by noting that the basic physics of transferring a narrow-bandwidth pulse 
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is very different from that involved in transferring an intra-cavity state. As the input 
state is centered around a single frequency, the transfer can essentially be viewed as a 
stationary scattering process from cavity- 1 input to cavity-2 output. As first noted in 
Ref. [15], high-fidelity transfer thus just becomes a set of requirements on the scattering 
matrix characterizing our system at a single frequency. The result is that high-fidelity 
can be achieved with a more modest set of parameters. One does not need to have 
many-photon strong coupling for each cavity (i.e. Gi ^ k)^ but (as we show) the more 
modest requirement of a large optomechanical co-operativity: Ci ^ 1 for each cavity, 
where 

a = ^. (32) 

In what follows, we solve for the scattering matrix describing our system, and use 
it to understand how and why one can achieve high-fidelity transfer of itinerant states. 
In particular, we demonstrate the mechanically dark-mode introduced in Sec. 2.2 also 
plays a crucial role in itinerant state transfer, allowing it to be effective in both regime 
of strong and weak optomechanical coupling. We consider both the transfer of Gaussian 
wavepackets, where again the simple form of Eq. (19) for the transfer fidelity holds, and 
the transfer of arbitrary non-classical states. 

4-1' Langevin-Heisenberg equation and scattering matrix 

The itinerant state transfer can be conveniently studied using the Heisenberg-Langevin 
equations of the whole system. Working as always in an interaction picture with respect 
to the main cavity drive frequencies, these take the form [32, 33]: 



a = - iujMd - - ^ X] Gidi - ^di, 



2 

di = iAidi - -^di - iGid - y^idi^.^^ - V^^^n (33) 

with di^i^ representing both input noise (taken to be white) and signals driving each 
cavity. In particular, the itinerant pulse incident on cavity 1 that we wish to transfer 
will be described by the operator di^in. K.i is the damping rate of cavity i arising from its 
coupling to the waveguide used to drive it and to extract signals from it. In contrast, is 
the damping rate of cavity i due to intrinsic losses, and d[ represents the corresponding 
noise. Finally, ain describes the thermal and quantum noise driving the mechanical 
resonator from its intrinsic dissipative bath. 

Solving Eq. (33) and using standard input-output relations [32] yield the relation 
between input and output fields 

lout[^] = s[uj]X^[uj] + [u] (34) 

with the vector of operators A defined as 

(^di[uj]j2[^id[u;]^ . (35) 
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1'= (d;M,4M). (36) 

Here, s is the 3 x 3[ujm] matrix describing scattering between the cavity 1 waveguide, 
the cavity 2 waveguide, and the thermal bath responsible for the mechanical dissipation. 
Note that in the general case where we have internal cavity losses, s will not be unitary. 
In contrast, s' is a 3 x 2 matrix which describes how noise from the sources of internal 
cavity loss can appear in the waveguides and in the mechanical bath. 

Both the above matrices can be found using Eq. (33) and the standard input-output 
relations 

AouM = A^^[uj] + (^, Vt) • A. (37) 
The elements are given by (z = 1, 2 denotes the two cavities) 

[uj] = 1 - 2/i, [uj] + aC [uj] /i,' [cj] . (38) 

Si2 [UJ] = 521 [^] = \/CiC2/il [Uj] /i2 [Uj] ( [^] , (39) 

i I — 

szi \oS\ = [cu] = - V Cii^i [u] C [oj] , (40) 

M = 1 - M , (41) 

with 

and fi [uj] is the susceptibility 



7 



^ -i{uj + Ai) ^-i{uj- ujm) 

and Ci the cooperativity of cavity i was defined in Eq.(32). 

The intrinsic noise scattering matrix s' [co] can also be calculated in the same way 




s.jMJ-^ (44) 



Assuming as always that large drive applied to each cavity (which yields the 
driven optomechanical couplings) are at the red-detuned sideband, i.e. = —^m^ 
the configuration of the scattering matrix is show in Fig. 9. As one can see, for both 
strong coupling (solid) and weak coupling (dashed), there is a transmission window 
peaked around frequency lum- Notice that the equations are in the rotating frame with 
respect to the drive frequency. Hence u — ljm corresponds in the lab frame to signals at 
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{oj-ojm)/G 

Figure 9. (Color online) The complex modulus of the scattering matrix |52i[co']| 
(Eq. (39) with no internal cavity losses, i.e. = 1) when cavities are in resonant with 
the mechanics = — cjm; S2i[uj] describes the transmission of a signal incident on 
cavity 1 to a signal leaving cavity 2 (and vice- versa). The two curves are plotted with 
the same mechanical damping rate 7. The solid/blue curve is obtained with hz/G = 0.1 
and G/7 = 10^. The red dashed curve is obtained with n/G = 2 and G/7 = 10^. 



the resonant frequency of each cavity. At this frequency the scattering elements become 

511 [uJm] =l-ri(^- ^^i^ ' (45) 

512 [^m] S21 [ujm] ^^^^^^ ^1^2 (46) 

ssi [^m] = Si3 [u)m] = ~ V ^, (47) 
2 

533 [u)m] = 1 - ^ (48) 



with 
and 



r^ = /^^/ (/^^ + Z^-) , (49) 



C = l + 4Ciri+4C2r2. (50) 



In the limit of symmetric cavities (i.e. = r and Ci = C2 = C)^ and high co- 
operativity C ^ 1^ the scattering matrix reduces to: 

r 

sn [u)m] = S22 [^m] ^l-r, 533 [ujm] = ^ ~ (^1) 
^512 [^m] = ^521 [ujm] ^ r (52) 

^531 [^m] ^513 [um] 523 [^m] ^32 [uJm] ^TJ^' (53) 

2v C 

We thus see that transmission at the cavity frequency becomes perfect in the limit of 
no internal loss (i.e. r ^ 1); in this limit, the above scattering matrix reduces to the 
results of [17]. In our following discussions of itinerant transfer, unless specified, we take 
this hmit for simphcity. 
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High fidelity transfer from di^in to rf2,out requires that over the input signal 
bandwidth, the transmission coefficient |52i[cc;]p ^ 1, as well as that |523[cl;]P ~ 
(i.e. negligible transmission of mechanical noise). Of course, if there is no internal loss, 
the scattering matrix s is unitary, and hence the former condition implies the latter 
condition. We see from the above results that in the limit of low internal cavity losses, 
these ideal scattering conditions are satisfied when the input signal is centered at ujm 
and the system parameters Ci — C ^ 1. In the following subsection, we give physical 
explanations of these perfect transmission conditions. 

4-2. Dark mode transfer due to Optomechanical EIT mechanism 

While the basic results for itinerant photon transfer follow immediately from the 
elementary calculation of the scattering matrix presented in the previous subsection, 
these calculations do not necessarily yield a good intuitive feel for why such ideal 
scattering conditions can be realized. Such a heuristic understanding is the goal of this 
subsection; in particular, we show that the mechanically-dark mode plays an essential 
role in achieving a high fidelity itinerant state transfer. 

To have protection against mechanical dissipation, one would ideally like the input 
state incident on cavity 1 to only excite the dark mode. Without dissipation, the dark 
mode Cdk defined in Eq. (9) is energetically separated from the mixed modes c± defined 
in Eq. (12). If the input signal is closely centered at frequency ujm in the displaced 
frame, only the dark mode is excited and hence protection against mechanical thermal 
noise is achieved. 

Including dissipation (and the consequent lifetime broadening), the input signal 
incident on cavity 1 will also excite the bright cavity mode Cbr (see Eq. (8)) as well as the 
mechanical mode a. This is not a issue in the strong coupling regime, as the transmission 
channel of the dark mode (central peak of solid line in Fig. 9) is well separated from 
the bright modes transmission (two side peaks of the solid line in Fig. 9). However, a 
problem arises in the weak coupling limit k> G\ the broadening of each mode is larger 
than their separation, and thus they are not resolvable from one another (see red dashed 
line in Fig. 9). Thus, for weak coupling, one cannot selectively drive only the dark mode 
by simply tuning the frequency of the input signal to be um- 

Fortunately, as we show below, in the weak-coupling regime, the unwanted 
excitation is irrelevant as long as the cooper at ivity of each cavity 1. In this limit, 

the bright mode amplitude (cbr) arising from the signal incident of cavity 1 is a factor 
^ l/C smaller than the dark mode amplitude, due to a destructive interference akin to 
the optomechanical analogue of electromagnetic-induced transparency [19, 20, 21]. 

In order to understand this, it is useful to introduce a "mixed" 3x3 susceptibility 
matrix x[^] which describes how the input modes incident on each resonator drive the 
non-local system modes C[uj\ — {cdkM, Cbr[ct;], a[cc;]}: 



C[uj] =x[a;]klinM, 



(54) 
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— * 

where k = diag{y^, V^} ^ diagonal matrix, and the vector is defined in 
Eq. (35). 

There are two main inputs driving this system: the signal from cavity 1 and the 
thermal noise from the mechanics. Hence the elements of particular interest are Xj,i 
Xj,a {j = {dk, br,a}). The matrix x[cl;] can be found directly from Eq. (33). While the 
full expression is rather lengthy, at frequency u — it can be written as (assuming 
no internal cavity losses) 



/ -G2M2 -GM -2i J'-"l. VC^2 \ 



XM = 



4Ci + 4C2 + 1 



-^2 
^2 

Gi G2 2i Gi-\-C2 



-2i^ -2i " 

7/^1 7/^2 



G2 1 



(55) 



with 

Consider first the excitation of the "bright" mode Cbr relative to the dark mode 
Cdk by a signal incident on cavity 1 at frequency uj — uom- In the ideal limit where 
Ci = C2 = C > 1, we find: 

Xbr,l [^m] 1 .g^x 

Xdk,i [^m] 4^ {ki + K2) 

Thus, the relative amplitude of the bright mode is suppressed by a large factor C 
compared to the dark mode, even if one does not have G ^ k and hence a strong 
spectral separation of these modes. Physically, this is due to a EIT-like destructive 
interference. In the absence of any optomechanical coupling (i.e. G = 0), an input 
signal incident on cavity 1 would excite both the dark mode and the bright mode equally. 
However, when G 7^ 0, the coupling between the bright mode and mechanics strongly 
modifies its susceptibility (this is the essence of optomechanical EIT). Heuristically, 
an excitation on the bright mode will be swapped multiple times between the bright 
mode and the mechanical mode, resulting in destructive interference. In case of two 
symmetrical cavities, the modified susceptibility due to these "multiple swaps" can be 
represented by a Dyson series: 

_ ^(0) ( 1 , ^2^(0) (0) , /^r2^(0) (0)\ , \ _ -^br,l _ Abr,l .^^x 

Xbr,l - Xbr,l \^ + ^ Xbr XV + [G Xbr X\ ^) + -J " ^ _ ^^2^(0)^(0) " ^ ^ 

where have focus on the ideal frequency u) = ujm- Here, XhJi ^ — \/2/a^i is the 
susceptibility of the bright (dark) mode to the input of cavity 1 when G ~ 0, 
and xi^^ = susceptibility of the mechanical mode, Xhr ^ —i/n\s the 

susceptibility of the bright mode to the mechanical mode in resonance. As we see, the 
coupling of the dark mode to the mechanical modes results in the suppression factor 
1/ (1 + G) in Eq. (58). This phenomenon identical to optomechanical EIT, where the 
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coupling of the optical mode with the mechanical mode prevents the cavity from being 
excited by an incident signal. The cancellation is valid as long as the the cooperativity 
is large, and independent of the strong coupling condition. 

On the other hand, the input signal on cavity 1 will also excite the mechanical 
mode a. Again using our mixed-susceptibility, one finds 

^^^-^./^^oc^ (59) 

Thus, if one is in the regime of weak coupling Gij^i < 1, the excitation of the 
mechanical mode can be appreciable. However, due to the relatively weak coupling 
7 between the mechanical resonator and its dissipative bath, the overall contribution 
from this excitation of a to Sout scales as y^a; this follows from the input-output relation 
^out = ^in + s/jo.. The net result is that only a small amount of the input signal is lost to 
the mechanical bath: the scattering matrix element 531 describing this process is small 
as \j\fC (c.f. Eq. (53)). Therefore, the transfer of the input signal thus occurs almost 
entirely via the dark mode in this limit. 

Good fidelity also requires that the dark mode, once excited by the input state, 
only leaks out via cavity 2, ensuring |52i[ct;M]| ~ 1. In another word, the input signal 
has little leakage through the output of cavity 1 . This requires a destructive interference 
between the promptly refiected input signal and the wave leaving the dark mode via 
cavity 1. That is, the signal refiected directly out of cavity 1 and the signal scattered 
into cavity 1 output cancels. For ^ 1, this interference cancellation results in the 
simple impedance matching condition [15, 17] Ci = C, i.e.: 

G\lKr = GIIk2 (60) 

We have thus shown how the incident signal primarily excites the mechanically-dark 
mode of the two-cavity optomechanical system, and how the dark mode then leaks out 
to cavity 2. It remains to show that we also have a suppression of noise originating from 
the mechanical bath in the output of cavity 2. Eq. (53) shows that for C 3> 1, ^23 \^m\ 
is suppressed by a small factor l/^fC] this directly results in the small value of the 
output heating parameter nh discussed in the main text (see Eq. (9)). This suppression 
is most easily understood in the case where a^i = /^2- In this case, the dark mode is 
completely decoupled from the mechanical mode. The bright mode is however driven 
by the mechanical noise, and thus provides a route for mechanical noise to corrupt the 
output from cavity 2. This process yields a contribution ^ v^'Xbr,a [^m] -^/t to 523 [cjm]- 
The relevant susceptibility Xbr,a [^m] ^ 1/G, yielding the result 523 [00 m] ^ \j\fC. 

The situation is only slightly more complicated when ki ^ In this case, the 
cavity decay terms in Eq. (33) can effectively mix the bright and the dark cavity modes. 
As a result, mechanical noise can first excite the bright mode, then be mixed into dark 
mode, then find its way into the cavity 2 output field. This additional process is not 
problematic, as it also scales the same way with C as the Ki = K2 process. This follows 
from the fact that Xdk,a[^M] ^ ((^1 - ^2)/^) Xbr,a[^M], where k = 2/^1/^2/(^1 + ^2) 
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(i.e. excitation of the dark mode by mechanical noise involves first exciting the bright 
mode). The net result is Eq. (53): 523[ct;M] scales as 1/VC irrespective of Ki/k2. 



4.3. Transfer fidelity for Gaussian states 

To quantify the transfer fidelity of the itinerant photon transfer scheme, we consider a 
Gaussian input state in a temporal mode (see, e.g., Ref. [33]) defined by 



A,in = (27r)- 



-1/2 



/ 



(61) 



where f[uj\ describes a wave packet incident on cavity 1 which is localized in both 
frequency and time; J rfcL;|/[cL;]p = 1 to ensure that Di^in is a canonical bosonic 
annihilation operator, di^i^ corresponds to an input mode with a fixed frequency, but 
completely delocalized in time. The fidelity of transferring this itinerant Gaussian 
state takes the same form as Eq. (19). In case of an itinerant coherent state input 
|'0in) oc exp {olD\ :^^ |0), the heating nh and the damping A can be explicitly written as 
(for details of the input temporal mode, see Appendix C) 

|2 



nh = ^ I du)\f[uj] S2i Ni 

2=1,2,M' 



A = \a\ min 



dwe-'^^S2i [u] 1/ [u] 



(62) 



We have optimized the final state pf in Eq. (18) over a time-translation r, so that if the 
output pulse is simply a time-delayed copy of the input pulse, F = 1. 

Taking our input mode to have mean frequency um and a Gaussian profile 

with variance Aoo^, and assuming Ci = C2 = C ^ 1, we find to leading order in Aoo: 



X 



M 



AC 



\a\ 



1 + 



~G 



1 f2Au\ 
8C ^ V~^/ 




(63) 



(64) 



Good fidelity requires a high cooperativity (7 Nm- In the weak-coupling regime 

G < one also needs Y^^lajAcc; <C (G^//^), which reflects the width of the S2i[uj\ 
transmission resonance. In the opposite regime G ^ Hi^ one needs Acj < K/y^|a| as 
shown in Fig. 10. 

In the limit Auj 0, Eq. (63,64) reduces to 

Nm 



AC' 



A ~ - — - 
8C 



(65) 



which shows the fidelity of itinerant state transfer depends solely on the cooperativity 
in case of zero-bandwidth. Further, we see that in comparison against the double-swap 
scheme, the mechanical-heating effect described by nh is reduced by a factor k/G. The 
expression of nh is the usual mechanical temperature expression of the cavity cooling in 
the weak coupling regime [23, 24]; unlike cavity-cooling, it describes nh in both weak 
and strong coupling regimes. 
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Figure 10. Fidelity F (red) and heating nh (blue) versus input bandwidth for 
transferring \a = Vs) coherent state. The red dash-dot line denotes the fidelity in the 
weak coupling case, where G/2tt = 0.1 MHz, k = 2tt x 0.5 KHz, and 7 = 27r x 10 Hz. 
The red solid line and the blue dashed line are both for strong coupling case where 
G/27r = 0.6 MHz, = 27r x 50 KHz and 7 = 27r x 1 KHz. The temperature is T = 2 K 
and the cavity frequencies are = 10 GHz, and /27r = 100 THz. 




4-4' Transfer of non- classical itinerant states 

Given the advantages of the itinerant transfer scheme, it is also interesting to consider 
how well it is able to transfer non-classical states. While in general it is difficult to 
obtain analytic expressions for the evolution of non-Gaussian states, we show that here, 
one can obtain useful and reliable analytic approximations. 

We again consider an input mode in a given temporal mode Di^i^] we take this 
mode to be centered on cjm, and for simplicity, to have a vanishingly small bandwidth 
Acj. Suppose now the input state incident on cavity 1 is 

Pm,l = 5Z<C^h)l,in(^l (66) 

mn 

where |n)i^in is a Fock state of this mode |n)i^in oc [D\ i^ |0). We also take the noise 
driving both cavities to be zero-temperature (A^i = A^2 = 0), but allow the mechanical 
resonator to be driven by thermal noise. Letting Pt\,{q^ A^m) be the probability of having 
q thermal quanta incident on the mechanical resonator, the fidelity can be decomposed 
as (see Appendix D for more details) 

00 q-r ^ 



g,r=0 mn d=—n 



with 



^ j\{j+d)\{n-j)\ir-j)\ Ws^s} ^' 
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Cooperativity C 



Figure 11. Fidelity versus cooperativity C = / in the narrow-bandwidth hmit. 
The input states are \a = ^/3) (green, dash-dot), |n = 3) Fock state (blue, dashed), 
\n = 1) Fock state (magenta, dashed) and (|1) + |3))/a/2 (red, solid). The coupling 
strength is G/27T = 0.6 MHz. Unless specified, the parameters are the same as the 
right column of Fig. 8. 



where 2F1 is the hyper-geometric function and s = s [ujm]- /n,o,t^'^~^~^ ^he amplitude 
for an input state |n, 0, q)^^ scattering into an output |r, n + rf, g — r — d)^^^. |n, 0, q) 
denotes the state with n photon in cavity 1, photon in cavity 2, and q phonon in the 
mechanics. Note that the scattering matrix conserves the total excitation, therefore the 
amplitudes between states with different excitation numbers vanish. 

In the case of a Fock state input |n), = Sjn,m the above result is reduced to 

00 00 00 

F = J2 Pir, n) = ^Y.P^i.{q, NJ {fj^^'^^f (69) 

r=0 r=0 q=0 

where P (r, n) is the probability of having r outgoing photons leaving cavity 1 and n 
photons leaving cavity 2, and 



(r,n) _ 




r 




S33 



511523 \ ^ / . ^ ^ . 531^23 \ 

2F1 [j -n,r-g;l+j; (70) 



.521513/ V 521533^ 

Note that in the regime of optimal state transfer Ci = C2 = C ^ 1^ the probability 
of having photons leave cavity 1 is small: the dark state effectively prevents mechanical 
photons from contributing, and Eq. (60) ensures minimal reflection of signal photons. 
One can thus get a good approximation by simply retaining the r = and r = 1 term 
in Eq. (69): F is approximately the probability of obtaining n photons in the cavity 2 
output mode and at most one photon leaving cavity 1. This is a rigorous lower bound 
on the exact fidelity, and is exact to order 1/C. 

In the limit C ^ 1^ one finds that to the leading order in 1/C the fidelity for 
transferring the n-photon itinerant Fock state is F 1 — [Nm (3 + 2n) + n] /AC. For 
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Nm ^ 1, the condition for a near-unity fidelity is thus C ^ NmTI] for a large-n Fock 
state, this is more stringent than the condition for having a large fidelity transfer of a 
coherent state with \a\ ^ yjn (c.f. Eqs. (63), (64)). 

The transfer fidelity of different non-Gaussian states together with a coherent state 
(for realistic parameters) are shown in Fig. 11. The figure shows high fidelity transfer 
of non-Gaussian states is possible if the cooperativity is sufficiently large. For the 
same cooperativity, the fidelity of transferring a Fock state (the blue dashed line) is 
lower than that of transferring a coherent state of the same average energy (the green 
dash-dot line). The difficulty probably arises from the non-trivial quantum nature and 
complicated phase space structure of the Fock state. 

5. Conclusion 

In this paper, we have studied quantum state transfer from an optical cavity to a 
microwave cavity via a jointly coupled mechanical resonator. We have introduced the 
notion of a mechanically-dark mode which is insensitive to the mechanical thermal noise. 
We have proposed different approaches to utilize the dark mode for state transfer and 
compared their advantages and disadvantages, in particular: 

1) For the transfer of an intra-cavity state, we find that a hybrid scheme which 
uses both the dark mode and the bright modes can achieve high transfer fidelity, even 
without pre-cooling the mechanical resonator to the ground state. We also provide a 
''phase diagram" to help experimentalists identify the best transfer scheme for their 
particular system. 

2) In the transfer of itinerant photon states, the dark mode can be used through 
a mechanism similar to the optomechanical electromagnetic transparency (EIT). 
Both Gaussian and non-Gaussian states can be transferred with high fidelity if the 
cooperativity is sufficiently high. 

In the short term, experiments are racing towards the transfer of coherent states and 
wave packets; we also expect the improvement due to dark modes can be demonstrated. 
Whereas in the long term, hopefully non-classical state transfer can be realized 
experimentally. This will open up new possibilities in hybrid quantum information 
network. 
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Appendix A. Gaussian state fidelity 

Since the state to be transferred is mostly a pure state, we consider the case where at 
least one of the states pi,pf is a pure state. Then the general definition of the transfer 
fidelity in Eq. (18) reduces to the simple overlap F = Tr (pipf). In terms of the Wigner 
representation W[ and Wf of each density matrix: 

F = Tr (piPf) = TT J^^ d^W, (e) W, (e) (A.l) 

Here, ^ is the vector formed by the quadratures of the modes of our system (A^ = 1 
here): <^ = {xi,pi,x,p2, • • • ,xn,Pn} and <^ = £2,^2, • • • .xn.Pn}- We define the 

quadratures of each mode as Xi — (^di + df^ /2 and pi = —i (^di — df^ /2, where di is the 
canonical annihilation operator of the mode. 

As our system is described by a quadratic Hamiltonian, an initial Gaussian state 
will remain Gaussian at all times. The Wigner function of a Gaussian state can be 
written in general in terms of the means and covariances of the mode quadratures: 



where V is the symmetrized covariance matrix, defined via Vjj' — | (^A^jAQf + AQfAQ)^ 

with Aij^ij-iij). 

Plugging Eq. (A. 2) into Eq. (A.l), we obtain 



F ^ 



(A.3) 
If we now define 

nh = 2Vdet(Vi + Vf) - 1 (A.4) 
.-((^>-(f,)).«^((^>-(f,)), (A., 

then Eq. (A.3) becomes the simple expression for the fidelity given in the main text, 
Eq. (19). Note that throughout the paper, we optimize the fidelity over simple rotations 
in phase space (so that if is a rotated version of pi, F — 1). 

Appendix B. Analytical solutions for intra-cavity state transfer 

For Gaussian states, we can rewrite the effective master equation in a general form 
characteristic to bilinear Hamiltonians 

-p (t) = -i [i- nl p {t)] +Y^^v (hli) p it) (B.l) 
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where H is a matrix form of the Hamiltonian, and the second term describes the effects 
of dissipation. We have three dissipative baths (one for each cavity, and the mechanical 
bath), and each can either absorb or emit energy, hence the index k runs from 1 to 6. 
For each term, describes the rate of either energy absorption or emission by one of 
the baths, and is the corresponding "jump" operator (which is hnear in the system 
quadrature operators). V denotes a standard Linblad superoperator: 

V (Lli) p{t)^2 iyj) p it) (Lli) ' - I (Lli) ' (Lli) , p (t) } (B.2) 

Using Eq. (B.l), one can get the time evolution of the average values (^) and the 
covariance matrix V: 



= Q(a^)>, -^r = QV + VQ^ + N (B.3) 



dt \ / dt 

with 

Q = 2<tH + 2(7 (Im r) , N = 2o- (ReF) (B.4) 
and the elements of a and F are defined as 



(Tij 



^mn — ^ Y^fc,m^M (B.5) 



Solving Eq. (B.3), the solutions are 



V (t) = e^' V (0) e^^' + f dr eQ(*--)NeQ^(*-") (B.6) 

Jo 

Using these results, one can calculate rzh and A in Eq. (A. 5) and the fidelity of state 
transfer. In the following, we show the result for double swap scheme, hybrid scheme 
and itinerant state transfer. As for the adiabatic passage scheme, high fidelity transfer 
requires a trade-off between the adiabaticity and the fast operation. The optimal 
evolution is thus not purely in the adiabatic limit; as a result, simple analytic expressions 
are difficult to obtain. 



Appendix B.l. Double swap scheme 



For the double swap scheme, the interaction is switched on sequentially. At each swap 
process, the relevant evolution and noise matrix can be derived from Eq. (B.4) to be 



— GiS 



^2 +7j_ _^ i^i 
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(B.8) 
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Here the subspace is the direct product of two systems under swap, i.e, in the first swap, 
the subspace is the direct product of the cavity 1 and the mechanical resonator; while 
the subspace of the second swap is the direct product of the mechanical resonator and 
the cavity 2. 

The initial state in cavity 1 to be transferred is assumed to be a squeezed state 



with 



\a,r)^D[a]S[r]\0) 



S [r] = exp ^dj - ^-df 



where r is taken to be real, and 

D [a] — exp ad\ — a'di 
Then the initial mean value and covariance matrix are 



Re a 
Im a 



1 



-2r 







,2r 



(B.9) 
(B.IO) 

(B.ll) 
(B.12) 



The initial state of the mechanical resonator is assumed to be pre-cooled into the ground 
state and the cavity 2 is also assumed to be in ground state initially. The swap time 
tis = ^/ (2Gz). Then the time evolution at = Yli'^i^ = 1^2) is given by Eq. (B.6) 
as: 

(Re a) e~ "^^^^^ 



6 its) = 



Vf (ts) = 



.-2r 





o2r 



(Im a) e ""'^'^ 



(B.13) 



4 \ e^^ / 4 

with Ki = {^i + 7) /4, z/^ = {ki - 7) / (4G^), A = Rijii + 70^1, B = ^2/^2 + 7/^2 and 

(ir (cos GiT + z/^ sin G^r)^ e~^*'^ 



(B.14) 



/ 



(ir (cos G^r — z/^ sin G^r)^e 



(B.15) 



where we have assumed the strong coupling limit and hence approximate 
— {n — 7)^ /4 G. Note that if pre-cooling is not performed, in Eq. (B.14) 
^1 ~^ ^1 (2A^M,o + 1)5 with A^M,o the initial phonon number of the mechanics. 

Plugging Eqs. (B.12), (B.13) and (B.14)) into Eqs. (A.3), (A.4) and (A.5), we obtain 
the full fidelity for transferring a coherent state. In the relevant strong-coupling limit 
where />:^,7 <C G, this results in Eqs. (23) and (24) in the main text. 

The calculation of the fidelity is only slightly more complicated in the case of a 
displaced, squeezed state, where now the parameter r > (c.f. Eq. (B.12)). Keeping 
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the phase of the displacement (i.e. a = lo^le^"^), the two parameters determining F 
become (keep until the leading order of (/^^ + 7) /G, also keep the higher-order terms of 
{hiiNi + ^Nm) /G as Nm can be large): 

nh[r] ^ ^1 + 2nh [0] cosh 2r + nl [0] + 2 (cosh 2r - l)^ Xiti^ - 1, (B.16) 



A[r] ^ A[0]^ 



cos^ 



-2r 



+ nh[0] 



+ 



sm 



.2r 



+ nh[0] 



(B.17) 



where nh[0], A[0] are heating and amplitude-decay parameters for a coherent state 
(c.f. Eq. (23)). Consider a squeezed ground state (i.e. Of = 0) in the strong squeezing 
limit r ^ 00. In this case, there is no amplitude-decay contribution to the fidelity 
A[r] = 0, and 

' (B.18) 



We see that the effects of mechanical heating are exponentially enhanced. Thus, for 
highly squeezed states (or other states with fine structures in phase space), a state 
transfer scheme which attempts to suppress these effects by using the mechanically- 
dark mode will be much better than the double-swap scheme. 



Appendix B.2. Hybrid scheme 

Both coupling Gi and G2 are switched on at the same time in hybrid scheme. The 
relevant evolution and noise matrix can be derived from Eq. (B.4) to be 



N= ^(1+Y) + ^(1-Y) 



(B.19) 



where we have taken Gi — G2 — G ^ tii — ti2 — i^-, 1 — 1 {2Nm + 1), ^ ^, 1 is a 6 x 6 
unity matrix and the definitions of S and Y are similar as Eq. (B.8) except that the 
dimension is increased to 6 x 6: 
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Then starting from a squeezed state as described by Eq. (B.12), the time evolution at 
ti,s = tt/v (with u = 2G^ — (/^ — 7)^ /16 ) is given by Eq. (B.6) as: 



_ /■v+7 



^hs 



( 



l + e 




(B.21) 
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V. (^.) = 4 — (l - ) + + - I; + X + 2 1 



with 



Jo Jo 

^ = / dre 2 sin z/r cos z/r, X— dre 2 '^cos z/r. (B.23) 
Jo Jo 

In the strong couphng hmit n <^G^ and considering a coherent state r = 0, the heating 

and damping can be simphfied to Eq. (31) in the main text. 

Appendix C. Transfer fidelity of itinerant Gaussian state 

In order to calculate the transfer fidelity for itinerant Gaussian state, we first define 
the input state to be generated by the mode creation operator An (An can be ^i,in, 
^2,in) which represents a time mode which spreads over a frequency interval Ao; around 
frequency um and localized in time around t = 

An = (27r)-^/' j dwf M 4 H (C.l) 

with rfin [uo] is the Fourier transform of the input operator which satisfies 
rfin [uJi\ , [uj2\ = 27rS (cji — 6^2), and f{uj) is normalized function which is localized in 
30th frequency and time. 

We define a smoothed operator [ujj] with discrete frequencies cuj = j • 6^ 

U e Z) [33] 

4 h] = dojd,, [^] ■ (C.2) 

Here 5^ is a small frequency interval on which / [uu] is approximately constant, and 
din i^j] satisfies the normal bosonic commutation relation 

d^ni^j] JL[^k]\ -Sjk^ (C.3) 
The integral of Eq. (C.l) is then changed into the summation 

An = E / 4 • 

j 

The output operator can be defined as 

Dout = (27r)-'/' J dujf [u] e-^('^-'^-)^(iout M (C.5) 
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with T is the time delay of the wave packet and a discretized summation can also be 
defined in a similar fashion as Eq. (C.4). 

Suppose the initial state is prepared in a coherent state 

p^\a)^{a\®\0)^{0\®Pa,th (C.6) 

with 

\a), = exp [aDl-^ - |0) . (C.7) 

The output state in Heisenberg picture is the same as Eq. (C.6), but the operators are 
changed from "in" to "out" according to the scattering matrix 

^^2,out M = S21 [UJ] rfl,in [^] + S22 [^] rf2,in [^] + "523 [^] ^in [^] (C.8) 

Using Eq. (C.l), (C.5), and taking the limit ^ and change the summation 
back to integral, we obtain the average value of quadrature in state Eq. (C.6) as: 

Re a \ / ^ \ j Re a( 



where 



Ci = (27T [ duje-'^^S2i M 1/ Ml') . (C.IO) 



And the covariance matrix are 

^"■'=i(i Vvu.= j(2»h + l)(j I I (C.ll) 

Using Eq. (A. 4) and (A. 5), we obtain the heating and amplitude damping parameters 
as shown Eq. (62) in the main text. 



Appendix D. Transfer fidelity of itinerant non-Gaussian state 

A Fock state input of cavity 1 is defined as 

Pl,in = ;^(^U)"|0)l(0|(l)l,in)" (D.l) 

where |0)^ is the vacuum state of the input mode and we only consider the input pulse 
with zero bandwidth for simplicity. We also assume the state of the cavity 2 is vacuum 
p2,in = 10)2 (0| , while the mechanical resonator is in a thermal state 

P<^,ir^-J2P'^^X^r.y\0)a{0\i^inr (D.2) 

q 

with pq =^ (1 — e~^^^) e~^^^. In terms of output operator, the output state is 

Pout = An = ^Pq \^q) {^q\ (D.3) 

q 
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with 



(D.4) 

and u = is the inverse of the scattering matrix. Hence the fidehty can be written as 



^yn\q\ 



q \ni,qa 



(D.5) 



According to Eq. (D.4), both un {o{l/C)) and u^i {o are small, i.e. all 

the processes to generate photon (reflecting) in cavity 1 are perturbations. One can 
therefore expand the state in terms of the photon numbers from the output of cavity 
1. Hence if we expand the fidelity as the summation of the probability of having r 
(r = 0, 1,2...) photon in the cavity 1 output, the series will converge quickly for large 
cooperativity. 

After some lengthy but simple algebraic calculations, we get the probability to have 
photon in cavity 1 output 

^13^32 \ 
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\2n I I 
Ul2\ \Uss\ 



2q 



Fi -n, 



U12U33J 



U33 



\2q-2 



Ul2\ 



\2n 



Ui2 ^ 



Fi (l-n,l-q;2;^^^ 
^ \ ' W12W33 



+K31- 2Fi(-n,i-g;i;tSt 



(D.6) 



(:d.7) 



where i^(2,i) is hypergeometric function. F^^^ is a lower bound of the exact fidelity, which 
physically corresponds to the case that the two cavities swap their states, so this fidelity 
represents a "swapping fidelity". Expansion of F^^^ in terms of s gives a contribution 
in the fidelity which is precise up to o (l/C^), 

F(^) ^ (1 - + g + ng) (D.8) 



The probability to have 1 photon in cavity 1 output is 



F33| 



U12 



\2n 



^u^. ,FWl-n,l-g;2;« 
+"31 • 2F, (-n, 1 - 5; 1; St) 



(D.9) 



^(0) _^ ^(1) jg ^ better lower bound of the actual fidelity. Physically it represents the 
fidelity of state transfer if the photon refiected from cavity 1 is measured to be or 1. 
Generally, to have r photons in the output of cavity 1, the corresponding probability 

is 
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The total fidelity is a summation of the all the probabilities F = F^'^\ In terms of 
scattering matrix, we get Eq. (69) in the main text. 

Furthermore, if the input state is a more general state pin^i — ^rnn^n^rn |^) (^L 
following a similar analysis, we can get Eq. (67) in the main text. 
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